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Abstract. We review and speculate on two recent developments of quantum optics and ultracold 
atoms. First, we discuss a possible realization of "phonon" physics, or as we call it refracton physics 
with optical lattices. To this aim we combine the physics of cold atoms with cavity QED, and 
investigate superfluid - Mott insulator quantum phase transition. The systems can exhibit cavity 
mode modifications due to local changes of refraction index (refractons). Second, we discuss the 
physics of strongly correlated particles in Abelian, and more interestingly in non-Abelian magnetic 
fields, using cold atoms. 



INTRODUCTION 

This paper has been presented as the final invited lecture at the International Conference 
on Atomic Physics, ICAP 2006, held in Innsbruck in July 2006. The speaker, M. 
Lewenstein has met all the advantages and disadvantages of being the last. The audience 
began to shrink and people were definitely starting to be tired of atomic physics, despite 
the enormous success of the ICAP 2006. On the other hand, being the last speaker M. 
Lewenstein had the opportunity to speculate in his lecture upon the future of atomic 
physics, and the physics of cold atoms in particular, and to thank truly the organisers. 

This paper consist of two relatively independent parts. In the first one we surmise 
the possibility of realising fascinating phenomena resembling those occurring in solid 
state physics due to the presence of phonons. We start by reviewing the superfluid-Mott 
insulator (SF-MI) transition of bosonic atoms in an optical lattice, which is a paradigm 
of a quantum phase transition. We turn then to the SF-MI transition inside an optical 
cavity. In the simplest situation, the atoms influence the cavity field by shifting the 
cavity resonance. The resulting effective Bose-Hubbard model suggest the existence of 
overlapping, competing Mott phases for wide range of parameters, and bistable behavior 
in the vicinity of the shifted cavity resonance. In a more sophisticated version, the model 
includes the possibility of local modifications of the cavity field due to the presence of 
atoms and local change of refractive index i.e. refracton physics. 

The second part deals with ultracold atomic gases in "artificial" non-Abelian fields. 
Again we review first distinct ways of realizing Abelian and non-Abelian magnetic field 
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in optical lattices using various types of laser manipulations etc. We present our recent 
results dealing with the problem of a single atom in a lattice, which, in the presence of 
Abelian (non-Abelian) magnetic fields, leads to a spectrum in the form of the famous 
Hofstadter butterfly (or yet not so famous moth). Finally, we speculate on non-Abelian 
fractional Hall effect with ultracold atoms. 



Introduction. Ultracold atomic gases in optical lattices (OL) are nowadays a sub- 
ject of very intensive studies, since they provide an unprecedented possibility to study 
numerous challenges of quantum many body physics (for reviews see [1, 2]). In particu- 
lar, such systems allow to realize various versions of Hubbard models [3], a paradigm of 
which is the Bose-Hubbard model [4]. This model exhibits supefluid (SF) - Mott insu- 
lator (MI) quantum phase transition [5], and its atomic realization has been proposed in 
the seminal Ref. [6], followed by the seminal experiments of Ref. [7]. Several aspect and 
modifications of SF - MI quantum phase transition, or better to say crossover [8], has 
been intensively studied recently (cf. [2, 9]). The control parameter of the SF - MI tran- 
sition is the strength of the optical potentials, i.e. laser intensity. For low laser intensities, 
bosonic atoms can easily tunnel around the lattice, and they condense into a maximally 
delocalized states with long range phase coherence, and density (on-site atom number) 
fluctuations. When laser intensity increases, optical potential barriers prevent tunneling. 
In effect, density fluctuations become costly, and the systems enters a gaped, incom- 
pressible MI phase, with no phase coherence, and reduced density fluctuations. 

The critique of pure perfectionism. Optical lattices provide an ideal toolkit for 
Hubbard models [3], since they are robust and do not exhibit phonons - in this sense 
they are perfect. As we often encounter in everyday life, perfectionism has positive 
and negative sides. In particular, the absence of phonons limits the whole richness of 
condensed matter phenomena that result from electron-phonon interactions, such as for 
instance the famous Peierls instability. Antiferromagnetic (non-frustrated) Heisenberg 
spin systems in perfect lattices, and in particular spin chains in ID at zero temperature 
often attain Neel ordering. In the presence of phonons, the system arranges in such a way 
that every second site move slightly to left, and every other second to the right, so that 
the energy of spin interactions is modulated: strong- weak-strong-.... The ground state 
has a form of a chain of dimers (singlets) located at the "strong" bonds. Contrary to Neel 
antiferromagnets which magnetize in arbitrarily small magnetic fields H, a dimerized 
antiferromagnet requires a a non-zero magnetization M to break the dimers. Antifer- 
romagnetic spin chains exhibit thus a magnetization plateau at M = 0. Magnetization 
plateaux with much complex origin have been recently intensively discussed in the lit- 
erature (cf. [10] and references therein). For instance, a simple XXZ spin chain, in an 
adiabatic approximation for "heavy" phonons is described by the Hamiltonian: 



TOWARD "REFRACTON" PHYSICS 
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Figure 1. Magnetization plateaux in a frustrated Heisenberg spin chain analogous to the one described 
by expression (1) (from [10]). 



where 5/ are lattice sites displacements, / and A nearest neigbour XY and Ising spin cou- 
plings respectively, and A describes the XY coupling modifications due to site displace- 
ments. This system is expected to exhibit novel magnetization plateaux at non-trivial 
rational values M = 1/3, 1/2.. (see Fig. 1). In the ground state, the displacement attains 
the self-consistent value 5 ; = (sjsi+i +sj +l st). It is challenging and fascinating to ask 
whether such physics can be captured and mimicked by cold atoms in lattices that are 
not perfect and react to their presence. A natural candidate for the realization of such 
systems is offered by Cavity Quantum Electrodynamics (CQED) [11]. 

Cavity QED and cold atoms. In fact, in CQED the atoms interact with the cavity 
mode and actively affect the cavity induced optical lattice. Several studies address the 
scaling of cavity QED dynamics with the number of atoms and their temperature. 
Cavity QED techniques were used to measure pair correlations in atom laser [12]. Self- 
organization of atoms in longitudinally and transversally pumped cavity was observed 
in [13], and theoretically characterized in [14]. Bragg scattering of atomic structures, 
which built up in the potential of the optical resonator, has been investigated in [15]. 
The interactions and cooling of a single atom in a quantum optical lattice in a cavity 
was studied in Refs. [16, 17]. Maschler and Ritsch [18] implemented the Bose-Hubbard 
model for atoms in the one-dimensional potential of an optical resonator, which was 
longitudinally pumped. They showed that the coefficients determining the dynamics, 
depend on the number of atoms and exhibit long range interactions, due to the photon- 
mediated forces. In a recent work, we have applied the model developed in [18] to 
study how the quantum phases of ultracold atomic gases are modified by the photon- 
mediated long-range interaction due to the optical resonator [19]. We focused onto the 
superfluid - Mott insulator phase transition for gases of few hundred atoms, and apply 
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the methods of quantum statistical physics, namely strong coupling expansion [20] to 
calculate the boundaries of the Mott insulating phases in the ji — r\ plane, where jl is the 
chemical potential, while r\ characterizes the pump strength. We predicted existence of 
overlapping, competing Mott phases, that may even consist of two not connected regions 
for a wide range of parameters. We predicted also dispersive bistable behavior [21] in 
the vicinity of the shifted cavity resonance in the strong coupling regime. 

Our model. Our model considers atoms confined in one dimension inside an optical 
resonator, which is pumped by a classical field. The atomic dipole transition is far-off 
resonance from the cavity mode, which induces a dipole potential on the atomic ground 
state. Using the notation of [18], the single-particle Hamiltonian reads: 

2 

Hq = ^ — h hUo cos 2 (kx)n^ — M c n p h — ihrj i^ci — . (2) 

Here, p and m are the atomic momentum and mass, r\ is the amplitude of the pump at 
frequency co p , A a = co p — CO a and A c = co p — (O c are the detunings of the pump from atom 
and cavity frequencies, g(x) = gocos(kx) is the atom-cavity mode coupling, k = (O c /c 
the mode wave vector, and a the creation and annihilation operators of a photon Ti(O c , 
«ph = c^a the number of photon, and Uq = gl/A a is the depth of the single-photon dipole 
potential. The many-body Hamiltonian is obtained from Eq. (2) in second quantization 
with the atomic field operators ^(x), and including the atomic contact interaction [18]. 
In the bad-cavity limit we eliminate, from the atomic dynamics, the cavity field variables 
and the number of photons takes the value 

k 2 + [Ac - U fdxcos 2 (kx)^(x)^(x)] 2 

with K the rate of cavity damping. The cavity potential depends thus non-linearly on 
the atomic density. In particular, certain atomic phases correspond to resonances, which 
increase the number of photons and thus the optical lattice depth. This is the salient 
physical property of this system, which gives rise to novel dynamics, as we will show. 
We now consider the regime of validity of the tight binding approximation (TBA) and 
expand the atomic field operators in the lowest energy band as ^(x) = Y*ibiw{x — xt), 
where bi is the atomic annihilation operator at site i, and w(x — xi) is the Wannier state 
localized around site i, which depends on the photon number n ph , and thus on the atomic 
density. In the Wannier expansion we keep on-site and nearest-neighbour couplings, and 
neglect all other couplings, and obtain the Hamiltonian 

H/U = -tB + \ Y^mim -\)-gB 2 - nN, (4) 

where N = Y*i n i = Y*ib]bi is the atom number operator, B = X^&J^-i +h.c. is the 
hopping term, and U denotes the strength of on-site interactions. The parameters 

f Ei | ifWi (* 2 -C 2 ) , rfhUffit (3k 2 - C 2 ) 

u u(k 2 +c 2 ) 2 ' 8 u(k 2 + ; 2 ) 3 



Travelling to exotic places with ultracold atoms 



5th February 2008 



4 



denote tunneling and long-range coupling respectively, and are expressed in terms 
of the Af-dependent term £ = A c — UoJqN and of the integrals Eg = f dxw(x — 
xi)(—h 2 /2m)(d 2 /dx 2 )w(x — xi + f) and 7<? = f dxw(x — xi) cos 2 (kx)w(x — xi + (>), with 
£ = 0, 1, where we have assumed 7i <C 7o. Note that in Eq. (4) the number of particles 
is conserved, since [N, H] = 0. 

SF - MI transition in a cavity. In [19] we have generalized and applied the method 
of H. Monien's group [20] (for more references see [2]) to calculate the lobes for a 
modified ID Bose-Hubbard model (4) describing atoms in an optical lattice created by 
pumping a laser beam into the cavity. Technically, this corresponds of using perturbation 
theory in B to calculate the ground state energy, and the energies of excited states with 
one additional atom (one hole) present and comparing them. The major difference to the 
standard case is that the intensity of the cavity field depends on the number of atoms 
present, since the atoms shift collectively the cavity resonance. Thus the coefficients 
t, U, g become very complicated functions of the cavity detuning, intensity of the 
pumping laser, N, etc. Moreover, quantum fluctuations of the resonance shift induce 
long range interactions between the atoms, proportional to g and hence the parameters 
t, U, g have to be calculated self-consistently. Moreover these calculations sometimes 
exhibit bistability effect! Instead of using /i and t to describe the phase diagram, it is 
necessary to use a proper control parameter which is the strength of the pumping laser. 
A typical phase diagram as a function of /i in the recoil units, and Jc/rj is shown in 
Fig. 2. The striking effect is the overlap a different Mott phases, and even a presence 
disconnected Mott regions, which follows from the fact that the expressions for t and U 
for n = 1, 2, 3, . . . Mott phases are different. 
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Figure 2. Overlapping and disconnected Mott insulator regions for A c = and Uq = 2k. 
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Finally, if A a and A c have the same sign, and one may encounter resonances such that 
A c — UoJqN = 0; the system may exhibit a bistable behaviour, as reported in [19] 

Toward genuine "refractons". An obvious challenging question is whether one 
can generalize the above results to a situation when the atoms affect the cavity mode 
locally. Such a regime of CQED is experimentally feasible and has been recently con- 
sidered by D. Meiser and P. Meystre [22], who termed it "the Super-Strong Coupling 
Regime of Cavity Quantum Electrodynamics". We present here a somewhat alternative 
approach based on the techniques that one of us developed together with R. Glauber, T. 
Mossberg and K. Zakrzewski to describe spontaneous emission and Lamb shift modifi- 
cations in dielectric media and/or in "colored" photonic reservoirs [23]. The idea is to 
use the full quantum electrodynamical description of the electric field inside the cav- 
ity, i.e. expand it into the continuum of modes that have a significant maximum in the 
density of modes at the cavity resonance. The single atom Hamiltonian becomes then 

H = ^ + W E\x)E(x)+ ! dkcka\k)a(k) 
2m J 

— ihrj J dkp(k) ^a(k)exp(— ic\k\t) — a* (k) exp (+ic\k\t)j cos (he), (6) 

where d(k), d'(k) are photon annihilation and creation operators of photons in standing 
wave modes (that inside the cavity may be well described by cos (foe)). The pumping 
laser pulse is expanded into these modes, and the amplitude of expansion, p(k) is 
strongly peaked at k ~ k^, and CO = c\k\ ~ c\ki\ = oo^. The electric field operator E(x) 
inside the cavity can be expanded as (for simplicity we consider only even modes): 

E( x ) oc / dk- ,a(*)cos(fcc) = qcos(*L*) 

J (k-co c + iK) K-i(A c -U 8(x)) 

where the final expression appears after adiabatic elimination of the cavity field. The AC 
Stark shift (which was given by 8 = JoN + JiB) becomes now local, i.e. ^-dependent: 

8(xi) = (exp(-K\xi - xj\) + exp(-K\xi+ Xj\) fij + J\B (8) 

j 

with hj being the atom number operator at the site j. 
In the hard boson limit we end up with the Hamiltonian: 

H= - J^[J(i,{ s z,j})siS i+ i+s] +l SiJ(i,{s z j}) 

i 

- E[a(/,{^})^+*M»"{^})] -E^, (9) 

Apart from the fact that the system is ferromagnetic, it resembles to great extend the 
spin chain with adiabatic phonons (1): the main difference being that the self-consistent 
hopping amplitudes depend functionally on the density (nfs, or equivalently s z j's), 

rather than on spin-spin correlations ((sjsj+i 
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ULTRACOLD GASES IN "ARTIFICIAL" NON-ABELIAN FIELDS 

Ultracod atoms and HEP. In the recent years the links and interconnections be- 
tween physics of ultracold atoms and condensed matter physics became solid and well 
established. Through these links, modern atomic physics reaches the frontiers of the 
quantum field theory, and has indirect relations to the modern high energy physics. Very 
recently, it has appeared several proposals on how to use ultracold atoms to simulate 
abelian U(l) lattice gauge theory [24]. Another research line investigates the possibility 
of inserting atoms in artificial non- Abelian magnetic fields by employing electromag- 
netically induced transparence [25], or optical lattices [26]. The recent progress in the 
latter case will be presented here. 

Artificial Abelian "magnetic fields" in a lattice. As it is well known, rapidly ro- 
tating harmonically trapped gases of neutral atoms exhibit effects analogous to charged 
particles in uniform magnetic fields (for a recent overview, see [27]). Thus one should 
be able to realize analogues of fractional quantum Hall effect (FQHE) in such systems. 
However, to achieve a regime of "strong" magnetic fields it is desirable to realize artifi- 
cial fields in an optical lattice by controlling e.g. the phases of the hopping amplitudes. 
In particular, single particle stationary Schrodinger equation in 2D square lattice in per- 
pendicular magnetic field, in Landau gauge (i.e. with the vector potential A = (0, Bx, 0)) 
reads: 

E^(m,n) = -J[ x ¥(m+l,n)+ x ¥(m-l,n)+exp(-i27tCCm) x ¥(m 1 n+l) 

+ exp(+i2jzccm) x ¥(m 1 n- 1)], (10) 

where (m, n) are integers enumerating coordinates on the lattice, and a — Ba 2 / 2% (he / e) , 
is the magnetic flux thought an elementary plaquette of area a 2 in units of quantum flux. 
Assuming plane waves in the y direction, ^(m, n) = exp(ivn)g(m) , we obtain the famous 
Harper's equation: 

£g(m) = - [g(m+ 1) + g (m - \)+ 2 cos(2n am -v)g(m)] . (11) 

Jaksch and Zoller [28] were the first to propose how to realize the model of Eq. (10) 
employing internal states of atoms, and using appropriate combination of laser assisted 
tunneling, lattice tilting, resonance pumping, etc. They argued that such systems could 
be used to study the Hofstadter butterfly (i.e. spectrum of states) in the presence of weak 
atom-atom interactions. Several other groups [29, 30] proposed since then alternative 
methods of realizing Abelian fields, aiming at the possibility of realizing FQHE in 
such systems: Laughlin states at low a [30], and analogues of bilayer FQHE (cf. [31]) 
Halperin states [32]. 

Artificial non-Abelian "magnetic fields" in a lattice. In [26] we demonstrated 
that using atoms with multiple internal states ("colours" or "lavours"), one can realize 
Non-Abelian magnetic fields in an lattice using a similar scheme as [28], i.e. employ- 
ing laser assisted tunneling, lattice tilting, resonance pumping etc. Generally speaking, 
the proposed scheme allows, among other things, to generate fields, whose vector po- 
tential is linear in x and y, i.e. is of the form A = (M x + N x (x/d) + O x (y/a),(M y + 
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Ny(x/a) + Oy(y/a),0), where M,N, and O are arbitrary hermitian matrices from the 
algebra of the considered group, which can be chosen practically at will: SU(n), 
U(n), GL(n) etc. For instance, with A = (m x ,Ny(x/a),0) with M x ,N y - hermitian, and 
N y = 2^diag( ai , (X2, • • • , OW) > we obtain the U (n) generalization of the Harper's equation 
for the rc-component wave function: 

£g(m) = — [exp(iM x )g(m+l)+exp(-iM x )g(m-l) + 2cos(2nN y m-v)g(m)] . (12) 

In Ref. [26] we have studied the spectrum of this equation for the SU (2) case, and 
rational d\ and «2- Bloch theory applies in such a case and the spectrum consists of 
allowed bands and gaps, similarly as in the Hofstadter problem for n = 1 and M x = 0. 
In the latter case the allowed energy band plotted against (rational) a, form the complex 
fractal figure known as Hofstadter butterfly. In the present case, we plot instead energy 
gaps. They form a fractal arrangements of genuine holes is the spectrum (Fig. 3), which 
we term Hostadter moth. 




Figure 3. Hofstadter moth's energy gaps plotted in various colours against rational values of "magnetic 
fluxes" (X\ and tti. Whole spectrum is contained between the two external, slightly curved, grey "walls". 

Despite its richness, Fig. 3 contain relatively little informations on density of states 
and corresponds to an infinite system. We consider realistic systems of 1000 sites size, 
and plot a section of the moth for a fixed value of, say, 0,2 = 2/5. Such a plot in which the 
red points (holes) indicate the presence (absence) of an energy eigenvalue is, perhaps, 
the simplest representation of the density of states (see Fig. 4). The moth section is 
especially impressive if one turns the figure 90 degrees clockwise. The moth looks like 
a horrible monster from a horrible Japanese cartoons, that are being shown to the kids 
all over the world in order to educate them in the spirit of war, killings, anger, etc. 

Non-Abelian fractional quantum Hall effect. We term by non-Abelian FQHE 
the analog of the "standard" (if one dares to use such adjective) FQHE, but in the 
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Figure 4. Section of the Hofstadter moth at ai = 2/5 in a lattice of 1000 sites. 



constant non-Abelian magnetic fields. Our starting point is a non-abelian SU(n) field 
in a 2D square lattice. Note, that we cannot use the gauge potential from the previ- 
ous subsection, since it leads to a non-zero gauge field tensor. Instead we use A = 
(—B(y — bSy)/2,B(x + bS y )/2,0), where S Xjy>z are the spin operators (S = (n — l)/2). 
This potential has a non-vanishing constant rotation, and is non-trivially non-Abelian, 
since its components do not commute. We consider now not too strong fields so that the 
continuum approximation in the lattice maybe performed. The single particle Hamilto- 
nian becomes thus for n = 2 



H 



(p x + eB(y-b6 y )/2c) 2 (p y -eB(x + b& x )/2c 



2m 



+ ■ 



2m 



+ coo z /2, 



(13) 



where we have included o z term to allow for non-trivial coupling of different Landau 
levels. Introducing the operators that annihilate or create Landau level number, D = 
(Q + iP)V2, D\ with Q = ^Jc/eB(p x + eBy), P = ^c/eB(p y - eBx), the Hamiltonian 
attains the form of the Jaynes-Cummings model [21]: 



H/hCl = IffD-ipid^D-ffd) 



i ~ £ T £ 

+ coo 1 o 



(14) 



where Q. = eB/hmc, (b = CO /CI, and /3 = b^/eB/2c. The many body Hamiltonian in- 
cludes the sum of single particle terms plus the sum of pair interactions of contact type, 



iPi&t&i-fij&i) +coala i + g -^8^(zr 

L ¥J 



(15) 



where n = Xj + iyt are complex coordinates in the plane. 

There are some similarities between (15) and a multicomponent quantum Hall system 
studied in the literature (see article by S.M. Girvin and A.H. Monald in [31]), such 
as electronic systems with spin, or bilayer systems. In the spin analogy our model 
introduces several novel aspects: i) we may use it for fermions or bosons, ii) we may 
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consider arbitrary spin (pseudo-spin), and in particular consider contact interactions of 
the spinor type (cf. [2]). Concerning multilayer analogy we stress that, similarly, we 
may consider more general "multilayer" =multicolour situations, and control coherent 
coupling between the "layers". Our model breaks the spin up-down symmetry, and 
necessarily mixes various Landau levels, if we allow for spin-up excitations. If all spins 
are down, the ground state (GS) takes the Lauglin form: 

\Q,E = 0} = Yl(z l -z J ) m cxp(-^\z l \ 2 /4)\ •-,!), (16) 

i<j i 

where 1/ra = v is the filling factor, E = £,Z)jZ) ; is the total number of excitations, and 
E its mean value. Correspondingly, the GS with one excitation will be 

\^\E=\)) = nfe-^rexp(-£|z ; | 2 /4)|W(T))|+£Z)Ja, J B=l) (17) 

i<j i i 

It is natural to generalize the concept of hole excitations for the non-Abelian FQHE, 
defining single hole states TLfe' — Z a ) | CO, E) , two hole states Y\i{u — Z a ){zi — Z b )\G),E) 
etc., where Z a , Z\, are n x n matrices. These excitations constitute obviously fractional 
anyons, but currently we are working on checking whether they themselves are non- 
Abelian, i.e. transform according the non-Abelian representations of the permutation 
group [33]). This has been a long lasting quest for non-Abelian anyons in the con- 
densed matter literature, but without no clear experimental observation so far. The most 
prominent candidate is electronic (fermionic) v = 5/2 FQHE state. This state has been 
observed in experiments, Moore and Read [34], and independently Greiter, Wen and 
Wilczek [35] proposed to explain it in terms of the "Pfaffian state" (see also [31]). Re- 
cently, however, Toke and Jain [36] proposed an alternative "composite fermions" model 
which does not relate to non-Abelian statistics in any obvious manner. For bosons, a 
promising candidate is the v = 3/2 state from the so called Read-Rezayi sequence of in- 
compressible correlated liquids. This state seem to be a true ground state of the rapidly 
rotating gas of bosons interacting via contact (Van der Waals) forces with a moderate 
amount of dipolar interactions [37]. Such situation may be achieved, for instance, with 
Bose condensed Chromium, as in experiments of the T. Pfau's group [38]. We hope 
that non-Abelian FQHE, due its profound and direct non-commutative character, will 
provide further, experimentally feasible examples of non-Abelian anyons. 

We conclude with the standard conclusion of M. Lewenstein and the motto to the 
review [2], which is a citation from William Shakespeare's "Hamlet": There are more 
things in heaven and earth, Horatio, than are dreamed of in your philosophy. In the 
present context, it expresses our neverending curiosity, enthusiasm, joy, and excitement 
of working in atomic physics in general, and physics the of ultracold atoms in particular. 
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port from German Government(SFB 407, SPP 1 1 16, GK 282, 436 POL), Spanish Gov- 
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G.M, C. M and J.L are supported by Spanish Ramon y Cajal Program, EC-Marie Curie 
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